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CHANDRA TEMPERATURE MAP OF ABELL 754 AND CONSTRAINTS ON THERMAL CONDUCTION 
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ABSTRACT 

We use Chandra data to derive a detailed gas temperature map of the nearby, hot, merging galaxy cluster 
A754. Combined with the X-ray and optical images, the map reveals a more complex merger geometry than 
previously thought, possibly involving more than two subclusters or a cool gas cloud sloshing independently 
from its former host subcluster. In the cluster central region, we detect spatial variations of the gas temperature 
on all linear scales, from 100 kpc (the map resolution) and up, which likely remain from a merger shock 
passage. These variations are used to derive an upper limit on effective thermal conductivity on a 100 kpc 
scale, which is at least an order of magnitude lower than the Spitzer value. This constraint pertains to the 
bulk of the intracluster gas, as compared to the previously reported estimates for cold fronts (which are rather 
peculiar sites). If the conductivity in a tangled magnetic field is at the recently predicted higher values (i.e., 
about 1/5 of the Spitzer value), the observed suppression can be achieved, for example, if the intracluster gas 
consists of magnetically isolated domains. 

Subject headings: Conduction — galaxies: clusters: individual (A754) — intergalactic medium — magnetic 
fields — X-rays: galaxies: clusters 



1. INTRODUCTION 

Since the non-detection by XMM and Chandra of the ex- 
pected quantities of cool gas in cluster central regions (Pe- 
terson et al. 2001 and later works), there has been a renewed 
interest in heat conduction as a possible mechanism to off-set 
radiative cooling (e.g., Narayan & Medvedev 2001; Gruzi- 
nov 2002; Fabian, Voigt, & Morris 2002). Detailed analyses 
of several cooling flow clusters (Voigt et al. 2002; Zakam- 
ska & Narayan 2003) show that heat inflow from the outer 
regions may indeed be sufficient over most of the cooling re- 
gion, but only if the conductivity, k, is not lower than Ks/3, 
where Ks is the Spitzer (1962) value. Such high conduction 
would have other important effects, e.g., erasing the cluster- 
scale temperature gradients (e.g., David, Hughes, & Tucker 
1992) and causing significant heat loss into the intercluster 
space (Loeb 2002). 

The conductivity of the cluster plasma should be reduced by 
the magnetic fields; however, theoretical estima tes for such 
a reduction have not yet converged (see §3 3.21 for details). 
On the observational side, David et al. (1992) derived a lower 
limit of k > Ks/10 for the Coma cluster. However, recent 
higher-resolution data revealed temperature variations in the 
cluster region that those authors assumed isothermal (e.g., 
Neumann et al. 2002), which may weaken their constraint. 
The only cluster-scale results derived so far using the high- 
resolution data were for the conduction across cold fronts. Et- 
tori & Fabian (2000) pointed out that the temperature jumps 
across the fronts in A2142 observed by Chandra (Markevitch 
et al. 2000) imply k < k.s/250. Although one can argue that 
as the cold front moves, the temperature gradient may be con- 
tinually sharpened because of stripping of the heated bound- 
ary layer, their conclusion should be correct. As shown by 
Vikhlinin, Markevitch & Murray (2001c) for A3667, there 
is no such boundary layer and the gas density discontinuity 
at the front is very sharp — much narrower than the elec- 
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tron mean free path, indicating that diffusion and conduction 
across the front are indeed suppressed. The likely reason is 
a layer of the magnetic field parallel to the front, which can 
form as the gas flows around the front and straightens the tan- 
gled field lines (Vikhlinin, Markevitch & Murray 2001b). An- 
other place where conductivity was estimated is the boundary 
between the hot cluster gas and two cool, galaxy-size dense 
gas clouds in Coma (Vikhlinin et al. 2001a). For these clouds 
to survive, conduction across their boundaries should be re- 
duced by a factor 30-100. 

Those measurements probe rather peculiar sites in clusters, 
i.e. boundaries between different gas phases which are ex- 
pected to have disjoint magnetic fields and so be thermally 
isolated. As we will show below, an estimate that is more 
representative of the bulk of the gas can be obtained from the 
temperature maps of hot merging clusters, such as A754. 

Previous optical and X-ray studies of A754 — a rich clus- 
ter at z = 0.054 — showed that it is the prototype of a major 
merger. It has a complex galaxy distribution, X-ray morphol- 
ogy and gas temperature structure (Fabricant et al. 1986; Bird 
1994; Slezak, Durret, & Gerbal 1994; Zabludoff & Zaritsky 
1995; Henry & Briel 1995; Henriksen & Markevitch 1996, 
hereafter HM96). It also exhibits a radio halo (Kassim et al. 
2001). In this paper, we use Chandra da ta to d erive a de- 
tailed gas temperature map for A754. In SI3 3.11 we dis cuss 
the merger details revealed by the new data. In 33 3.21 our 
temperature map is used for an estimate of thermal conduc- 
tivity of the cluster gas. At this redshift, 1 " = 1.13 kpc for our 
assumed Hq = 65 h^ km s" 1 1 



'Mpc" 



2. DATA ANALYSIS 

A754 was observed by Chandra in October 1999 using 
ACIS-I 3 for a clean exposure of 39 ks. Our analysis proce- 
dure, briefly outlined below, follows that in Markevitch et al. 
(2000) and Markevitch & Vikhlinin (2001) with a few up- 
dates. For background subtraction, we used the ACIS blank- 



3 Chandra Observatory Guide http://asc.harvard.edu/udocs/docs/ "Obser- 
vatory Guide," "ACIS" 



2 




137.4° 137.3° 137.2° 5 10 

« T, keV 



FIG. 1 . — Projected temperature map of A754 (colors) overlaid on the ACIS 0.8-5 keV image (contours, spaced by V2). Point sources are excluded. The map 
is adaptively smoothed (as part of the fitting procedure) by a Gaussian with a between 25" and 55"; the image is smoothed with a between 10" and 14". In 
most areas, different colors correspond to significantly different temperatures. Direct temperature fits for regions shown by white rectangles are given in the right 
panel (errors are 90%). The vertical dotted line shows the cluster average temperature. The step of the rectangles is 100 kpc. 



sky dataset normalized by the ratio of the 10-12 keV rates 
in the data and the background files (see Markevitch et al. 
2003); this normalization was within 1% of their exposure ra- 
tio. From the ROSAT All-Sky Survey (Snowden et al. 1997), 
there is no soft Galactic excess at the cluster position so no 
additional background modeling was necessary. The detec- 
tor low energy response was corrected for the secular change 
of the ACIS quantum efficiency (Plucinsky et al. 2002). All 
point sources were masked out. 

A single-temperature fit to the spectrum for the whole clus- 
ter (as far as we could reach, r = 9' = 0.6 Mpc) gives 
T e = 10.0 ±0.3 keV and a metal abundance of 0.30 ±0.05 
(relative to Anders & Grevesse 1989) at the 90% confidence. 
This temperature is somewhat higher than the ASCA value of 
8.5-9.0 keV (HM96); however, ASCA covered more of this 
highly nonisothermal cluster. The absorption column, if fitted 
as a free parameter, is in good agreement with the Galactic 
value of Nh = 4.36 x 10 20 cm" 2 , and therefore was fixed. 

To derive a temperature map, we extracted cluster images in 
7 energy bands between 0.8-9.0 keV, excluding the intervals 
1.8-2.2 keV because of the poor calibration and 7.4-7.6 
keV containing a bright fluorescent background line. These 
images were adaptively smoothed (identically in all energy 
bands). In each image pixel, the 7 flux values were fitted by 
a thermal model with fixed Nh and abundances. The result is 
shown in Fig.[Q 

3. DISCUSSION 

The large-scale temperature structure is in broad agreement 
with earlier ROSAT and ASCA results (Henry & Briel 1995; 
HM96). The new map offers a much more detailed look into 
this interesting cluster. Of note are a large hot area to the south 



and southwest (already seen in the earlier data), cool gas at the 
northeast tip of the elongated central bright body, and numer- 
ous small-scale temperature variations. Most of these varia- 
tions are statistically significant; to illustrate this, we extracted 
and fitted the real spectra for several less-convincing features 
(right panel in Fig.[0. There is also a curious wavy structure 
in both the X-ray brightness and the temperature map (in the 
general area of regions 1-6 in Fig. 0, whose origin we find 
hard to explain. 

3.1 . Merger geometry 

HM96 proposed that A754 is an off-center merger, and 
Roettiger, Stone, & Mushotzky (1998) succeeded in reproduc- 
ing the general large-scale X-ray structure in hydrodynamic 
simulations. In Fig. 2, which overlays X-ray and optical im- 
ages, one can discern two main galaxy concentrations, one 
around the western brightest galaxy, another just south of the 
X-ray brightness peak (e.g., Fabricant et al. 1986; Zabludoff 
& Zaritsky 1995). While the western X-ray extension (not 
completely covered by our data, but see, e.g., Henry & Briel 
1995) is associated with the western galaxy concentration, the 
second galaxy subcluster is offset from its presumed X-ray 
counterpart, the bright elongated core (Zabludoff & Zaritsky 
1995). In the Roettiger et al. simulations, this elongation and 
offset are caused by the ram pressure exerted from the south- 
east by the gas that used to belong to the western galaxy sub- 
cluster (which is now moving further to the northwest). 

However, the Chandra data complicate this picture by re- 
vealing interesting structure within the elongated gas core. 
There is a steep X-ray brightness drop-off at the northern 
tip of the body (Figs. 1, 2) which, although not as sharp 
and regular-shaped as some cold fronts (e.g., Vikhlinin et al. 



TEMPERATURE MAP OF A754 AND THERMAL CONDUCTION 



3 




FIG. 2. — Chandra X-ray brightness contours (same as in Fig. 1) overlaid 
on the Palomar Digitized Sky Survey optical image. 



2001b), indicates a similar moving interface. The elongated 
cloud apparently moves to the north (in projection) from the 
big galaxy concentration, not being swept away by ram pres- 
sure from the south, but as if dragged along by another sub- 
cluster. The cold region at the nose of the elongated body ap- 
pears to be the former core of that subcluster. However, it does 
not coincide with any large galaxy concentration (although 
there is a possible small group there, which alone would not 
be considered statistically significant). Alternatively, in the 
course of the merger, the elongated cloud may have com- 
pletely decoupled from its former dark matter host and is now 
"sloshing" independently. Such a state would be very short- 
lived, because a moving gas cloud without the dark matter 
gravity support is quickly destroyed by gas-dynamic instabil- 
ities. 

The coolest spot does not coincide with the gas density 
peak, so apparently it is not in pressure equilibrium with the 
rest of the elongated core. Such a transient state may arise, 
for example, if the cool core travels down the strong ambi- 
ent pressure gradient faster than its internal sound speed. Its 
northern tip may then expand too fast for the internal pressure 
to equalize. It is also possible that projection effects are at 
play (e.g., the elongated core is a projection of two clouds). 

We finally mention a brightness feature that looks like a 
bow shock leaving the merger site, first noted in the ROSAT 
image (Krivonos et al. 2003). It is located at the leftmost 
brightness contour in our figures (unfortunately, we cannot 
study it because of the limited coverage). It is consistent with 
the general merger picture proposed by Roettiger et al. 

To conclude, simple scenarios cannot explain all the com- 
plex details in the new X-ray data on A754. However, the 
very complexity of the temperature map of this cluster can be 
used for an interesting physical estimate, as we show below. 

3.2. Limit on thermal conduction 

In a plasma with a temperature gradient VT e such that the 
gradient's linear scale lj = T e /\WT e \ X e (where X e is elec- 



tron mean free path), heat flux is given by q = —K\7T e . In 
the absence of magnetic fields and for typical cluster plasma 
parameters, the collisional conductivity is given by Spitzer 
(1962): 

k s = l.3ln e X e k(kT e /m e ) l/2 (1) 

where In A is the Coulomb logarithm. For magnetic fields ex- 
pected in clusters (B > 0.1/iG, Carilli & Taylor 2002), the 
electron and proton gyro-radii are many orders of magnitude 
smaller than their Coulomb mean free paths. Therefore, con- 
duction occurs predominantly along the magnetic field lines, 
and the effective conductivity is determined by the topology 
of the field. If the field is chaotic with a coherence scale 
h 3> X e , the conductivity on scales lj ^ lg is reduced by a 
factor of 3 (e.g., Sarazin 1988; Gruzinov 2002). If l B < X e 
(which is likely the case — radio data suggest lg ^5-30 
kpc, e.g., Carilli & Taylor 2002, while X e ~ 10- 15 kpc in 
the central region of A754), the effective conductivity is re- 
duced because electrons have to travel a longer path along 
the tangled field lines to diffuse over a given distance. The 
reduction is of order (It/Ib) 2 (e-g-, Tribble 1989; note that 
k becomes scale-dependent), but probably not more than a 
factor of 100 if one considers the physical origin of the tan- 
gled field (Rosner & Tucker 1989; Tao 1995). In addition, the 
conductivity along the field lines may be reduced by a fac- 
tor of ~ 10 (for Ib ~ X e ) by magnetic mirrors in the stochas- 
tic field (Malyshkin & Kulsrud 2001). On the other hand, 
for certain spectra of the field fluctuations, diffusion across 
the field lines and divergence of the lines may boost the ef- 
fective conductivity to within a factor of only ~ 5 below ks 
(Narayan & Medvedev 2001). This requires a field tangled on 
a range of scales down to ~ lg / 100. We note that it is unclear 
whether such a field would be consistent with the observa- 
tions of the ~ 10 kpc coherent "patches" of polarization of 
the background radio sources (Carilli & Taylor 2002; Govoni 
2001), if those measurements indeed probe typical intraclus- 
ter fields. To summarize, the collisional thermal conductivity 
on 100 kpc scales is expected to be reduced by a factor 3- 100. 

Our temperature map of A754 offers an interesting oppor- 
tunity to put observational constraints on n. All over the 
cluster, the temperature varies on linear scales starting from 
100 kpc (the map resolution, limited by statistics) and up. 
This scale is so small that the temperature gradients should be 
erased very quickly if the conductivity is at the Spitzer value 
0. An approximate expression for the conduction timescale 
fcond = -(dlnT e /dt)~ l can be found in, e.g., Sarazin (1988): 

fcond ~ kn e lj / k s ^ 1.2 x 10 7 yr 

n. ( _h_ Y( kT * V 5/2 (2) 
X 2 x 10" 3 cm" 3 V 100 kpc ) \ 10 keV J ' 

-3/2 

where n e is from Jones & Forman (1999); f con d oc h 65 . Note 
that lj as defined above and determined from the map in Fig. 
1 would be greater than 100 kpc. However, we measure tem- 
peratures in projection, with a large column of gas at the mean 
cluster temperature on the line of sight, so the true tempera- 
ture variations should be significantly stronger. 4 We therefore 

4 The average true distance between regions of different temperature in the 
map is also greater, but it is the size of each region, not the distance between 
them, that determines It . 
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chose to use the apparent scale of the variations in the map 
as a realistic value of It for this estimate; this, along with the 
strong dependence of Ks on T e , are our greatest sources of un- 
certainty. We did not take into account mild saturation of the 
classical heat flux (for our X e and It one expects k ~ 0.7 ks, 
cf. Cowie & McKee 1977 and Sarazin 1988), which in reality 
should not occur because of the tangled fields. In any case, it 
would not affect the conclusion. 

At the same time, we do see these temperature variations 
everywhere in the cluster, so they should have lived for at least 
as long as it takes for a merger shock and/or a disturbing sub- 
cluster (that presumably have created them via compression 
and gas-dynamic instabilities) to travel across the central re- 
gion of the cluster, L > 800 kpc (approximately the size of 
our map). The sound speed in this region is a safe estimate of 
that shock's velocity regardless of its Mach number, since this 
region now contains the post-shock gas. This timescale is 



^age 



L 



5 x 10 s 



800 kpc VIOkeV 



kT e 



-1/2 



yr; (3) 



f age oc /igj. Comparison of eqs. (|2j and (|3} shows that the con- 
ductivity should be reduced by a factor 



(k/k S ) 1 ~ fageAcond > 40 ft, 



1/2 
65 ' 



(4) 



although we emphasize that, given all the uncertainties, this is 
only an order-of-magnitude estimate. 

The limit in eq. @ is inconsistent with the upper end of 
the theoretical range for k. We note that the "local" conduc- 
tivity considered by theorists can still be higher, if, for ex- 
ample, the ICM consists of thermally isolated domains with 
self-contained field structure. Such domains of gas that be- 
longed to different subclusters may survive if mixing during 
the merger is not very efficient, as is indeed suggested by the 
existence of cold fronts (Vikhlinin et al. 2001b; Narayan & 
Medvedev2001). 

Our upper limit on conduction appears to exclude it as a 
mechanism to off-set cooling in the cluster cores (where one 



needs n > Ks/3, see refs. in §1). However, conduction in a 
denser, cooler core is likely to operate in a different regime; 
for example, there one expects X e <C lg compared to X e ~ lg 
in the regions that we studied above. In addition, the central 
magnetic fields may have a different topology because of the 
competing effects of the radial cooling inflow (Bregman & 
David 1988; Soker & Sarazin 1990) and radial infall of the 
subcluster pieces vs. the widespread gas "sloshing" (Marke- 
vitch, Vikhlinin, & Forman 2002). This may result in a dif- 
ferent reduction factor. Finally, we note that our limit in eq. 
@, as well as the cooling flow conduction estimates, are rel- 
ative to Ks which itself is much lower in the cool cores. Our 
absolute upper limit on conductivity is not necessarily lower 
than the absolute conductivity needed to compensate cooling. 
Thus, even if the reduction factor for the Spitzer conduction is 
the same everywhere in the cluster, we cannot exclude a suf- 
ficiently high heat influx into the cooling core that is due to a 
different physical process (and so has a different dependence 
on gas parameters), such as turbulence, for example. 

4. SUMMARY 

We have derived a detailed temperature map of the proto- 
typical merging cluster A754. Simple merger scenarios can- 
not reproduce the complex details revealed by the new data; 
either A754 is a three-body merger, or the cool dense gas 
found in the cluster's bright elongated core has decoupled 
from its former subcluster host and is sloshing independently. 
We use the small-scale structure in the temperature map to 
derive the first constraint on thermal conductivity in the bulk 
of the gas in a cluster. The conductivity on scales ^100 kpc 
appears reduced by at least an order of magnitude from the 
Spitzer value. 
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